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REALIZATION PROBLEMS FOR LIMIT CYCLES OF PLANAR 
POLYNOMIAL VECTOR FIELDS 

JUAN MARGALEF-BENTABOL AND DANIEL PERALTA-SALAS 


Abstract. We show that for any finite configuration of closed curves F C 
one can construct an explicit planar polynomial vector field that realizes F, 
up to homeomorphism, as the set of its limit cycles with prescribed periods, 
multiplicities and stabilities. The only obstruction given on this data is the 
obvious compatibility relation between the stabilities and the parity of the 
multiplicities. The constructed vector fields are Darboux integrable and admit 
a polynomial inverse integrating factor. 


1. Introduction and statement of the main theorem 

We consider the planar vector field 
(1.1) X = P{x,y)d^ + Q{x,y)dy , 

where P and Q are polynomials. The degree of X is defined as the maximum 
of the degrees of P and Q, and will be denoted by deg(A). The study of the 
limit cycles of planar polynomial vector fields has a long tradition, starting with 
the celebrated second part of Hilbert’s 16th problem, cf. [12], which consists in 
finding the maximum number of limit cycles for the vector field O) in terms of 
the degree deg(A), and studying the relative positions of these cycles. We recall 
that a limit cycle of A is a periodic trajectory that is isolated. Despite having been 
formulated more than a century ago, Hilbert’s 16th problem is still wide open, even 
for quadratic vector fields, i.e. deg(A) = 2. 

A related problem that has attracted considerable attention in the last years is 
the realization problem for limit cycles. To state it in a precise way let us introduce 
some basic definitions. Let C be a closed curve embedded in A configuration 
of cycles is a finite set F = {Ci,..., Cn} of closed curves. 

Definition 1.1. We say that two configurations of cycles r,r' are equivalent if 
there exists a homeomorphism </> : —>■ such that (/>{r) = F'. A vector field X 
is said to realize a configuration of cycles F if its set of limit cycles is equivalent to 

F. 

The realization (or inverse) problem asks if, for any configuration of cycles F, 
there exists a planar vector field X that realizes F. In this problem it is usual to 
prescribe other dynamical properties of the limit cycles (e.g. stability) as well as 
some additional conditions on the vector field X (e.g. regularity). 

The realization problem for limit cycles was first addressed by Al’mukhamedov [1] 
for vector fields using the theory of Lyapunov functions. The same techniques 
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allowed Sverdlove [18] to solve the realization problem completely (also pre¬ 
scribing the stability of the cycles), fixing a gap in APmukhamedov’s construction 
that had been noticed by other authors, see [TH] for details. 

Regarding polynomial vector fields, the first realization result was obtained by 
Bautin in |3] (later corrected in 0 ), where explicit expressions for polynomial 
planar vector fields with a prescribed set of algebraic limit cycles were derived. 
Alternative constructions appear in the works of Winkel m, Christopher [^ and 
Korchagin M- The idea of using inverse integrating factors to construct polyno¬ 
mial vector fields with prescribed algebraic limit cycles was first given in [nmi. 
This idea culminated in the complete solution to the realization problem in the 
polynomial setting obtained by Llibre and Rodriguez in [15j using inverse integrat¬ 
ing factors and the Darboux theory of integrability. An alternative proof using 
Lyapunov functions was given in with a generalization to higher dimensions. 
A remarkable recent solution of the realization problem was obtained by Coll, Du- 
mortier and Prohens using polynomial Lienard equations and the theory of 
slow-fast systems. All these works provide an upper bound for deg(A) which de¬ 
pends on the configuration of cycles. Notice that since the degree of X is not fixed 
a priori, the realization problem is much easier than Hilbert’s 16th problem. 

The main result of this paper is a realization theorem for planar polynomial 
vector fields prescribing not only the configuration of limit cycles, but also their 
periods, multiplicities and stabilities. These are the three most basic invariants 
under smooth conjugacy, so it is reasonable to consider them as a part of the 
realization problem. The limit cycles of the systems constructed in m and |5] are 
all hyperbolic and their stabilities depend on the configuration that is realized (of 
course, there are no semistable limit cycles due to the hiperbolicity). 

To state our main theorem let us introduce some notation. We call Dc the 
compact set bounded in by the closed curve C. We say that a limit cycle C is 
stable (unstable) in the interior if all the trajectories of X in Dc that are sufhciently 
close to C approach the limit cycle as t —>■ oo (f —>■ —oo). Analogously, if all the 
trajectories of X in \ Dc that are close enough to C approach the limit cycle 
as t —>■ oo (t —>■ —oo), we say that C is stable (unstable) in the exterior. It is 
obvious that the exterior stability of a cycle is determined by its interior stability 
and multiplicity. 


Theorem 1.2. Consider sets {Ti,... ,T„} of positive constants and {mi ,..., m^} 
of positive integers. Let T = {Ci,..., C„} be a configuration of cycles with fixed 
interior stabilities. Then T is realized as the set of limit cycles of a polynomial vector 
field X, where each Ck has multiplicity m^, period and the required stability. 
Moreover, there exists an explicit upper bound for the degree of X in terms of the 
prescribed quantities (see Theorem 5.2). 


The proof of Theorem |1.2| is based on the construction of Llibre and Rodriguez m- 
As in m, the realized limit cycles are algebraic and the vector fields are Darboux 
integrable and admit a polynomial inverse integrating factor. We would like to re¬ 
mark that the limit cycles of a polynomial vector field do not need to be algebraic, 
cf. |16j . but we are not aware of any construction of polynomial vector fields with 
given limit cycles that are not algebraic. It thus remains an interesting open prob¬ 
lem to characterize those analytic curves that can be limit cycles of a polynomial 
vector field. 
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The structure of the paper is as follows. In Section we construct a vector field 
Xt realizing the configuration of cycles T with prescribed periods. The construction 
of a vector field Xm with prescribed limit cycles and multiplicities is presented in 
Section where we also combine both constructions to obtain a vector field Xxm 
prescribing periods and multiplicities at the same time. In Section we construct 
a vector field Xts realizing the configuration of cycles T with prescribed periods 
and stabilities. Finally, using the constructions in the previous sections, the main 
theorem is proved in Section 

2. A PLANAR POLYNOMIAL VECTOR FIELD WITH PRESCRIBED LIMIT CYCLES 

AND PERIODS 

In the following theorem we show that any configuration of cycles T can be 
realized by a planar polynomial vector field, where the period of each limit cycle 
is also prescribed. The method of the proof is a variation of the construction 
introduced in US]. In the statement of the theorem, a curve Cfc S F is called 
primary if no other curve Cj S F, j k, is contained in the domain Dck ■ 

We recall that a Darboux first integral is a (possibly multivalued) function G 
that can be written as: 

L 

1=1 

where fi, g and h are complex polynomials and A; G C are complex constants, and 
a smooth function V is an inverse integrating factor of a vector field X if 

X -XV = VAiY{X), 

where div denotes the divergence operator. 

Theorem 2.1. Let F = {Gi ,..., Gn} be a configuration of cycles and {Ti,..., T^} 
a set of positive constants, then F is realized by a planar polynomial vector field 
Xt with deg(AT) < 2(n + r), where each periodic orbit Ck has period Tk- Here r 
is the number of primary cycles in F. Moreover, Xt admits a polynomial inverse 
integrating factor, it is Darboux integrable and all its limit cycles are algebraic and 
hyperbolic. 

Proof. It is well known m that any configuration of cycles is homeomorphic 
to a set of circles. Accordingly, we can take that F consists of n disjoint circles of 
radii {rk}k^i centered at the points {pk}k=n i-®- 

(2.1) Gk = {fk{x,y) = 0} with fk{x,y):={x-Xkf + {y-yky-rl, 

where pk = {xk,yk)- We can safely assume that the first r cycles {Ci,..., Cr} are 
primary and that no point pk is contained in any circle Gj . 

Let us introduce the following auxiliary functions: 

9 k{x, y) := (x - XkY + {y - ykY , 

n 

A:=Y[fk, 

k^l 
n 

At :=n/"’ 

/c=l 


:= n gk , 




C := exp ( -2^6 »a: ) , 


k^l 


Dj' := A.j'BC , 
H'j' i = In Dy , 
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where are positive constants that will be fixed later in order to prescribe the 
desired periods Tj,, and 9^. is an angular (multivalued) function defined as 


( 2 . 2 ) 


9k ■= arctan 


y-Vk 
X- Xk 


We claim that the vector field Xt in the statement of the theorem can be defined 
as 

DHt 


Pt := -AB 


Xt '■= Prdx + Qrdy 


with 


Qt := AB 


dy 
BHt 
dx 


We observe that in the case that = 1 for all fc, then At = A and the vector field 
Xt is the same as in the construction by Llibre and Rodriguez m- Throughout 
the paper, we shall use the notation Xtr ■= XT\Tk=i- 

To prove that Xt satisfies all the claims in the statement of the theorem, let us 
first show that it is a polynomial vector field with deg)^^-) < 2(n + r), that Dt is a 
Darboux first integral and that Vt '■= AB is an inverse integrating factor. Indeed, 
noticing that 


(2.3) 


dx dy ’ 


^d_C^ d_B 

dy dx 


a straightforward computation shows that Pt and Qt can be written as 


(2.4) 

(2.5) 
where 




k^l 


^ik ■■= n ■ 

j^k 

It is clear from these expressions that Pt and Qt are polynomials of degree at most 
2(n + r) — 1. Moreover, it is easy to check that the Darboux function Dt satisfies 
Xt ■ XDt = 0, thus implying that Dt is a Darboux first integral of Xt- Another 
easy computation shows that the vector field ^ is divergence-free in \ 14^7 ^(0), 
and hence Vt is an inverse integrating factor of Xt- 

Accordingly, all the limit cycles of Xt are contained in the zero set Vijr^(O) of 
its inverse integrating factor m- Therefore, since 

R-i(O) =ru{pi,...,pr}, 

we conclude that if Xt has a limit cycle, it has to be precisely one of the circles 
{Cl,..., C„}. Let us now prove that indeed all of them are realized as limit cycles. 

We first show that each Ck is a periodic trajectory of Xt- Since Vljr^(O) is 
invariant under the flow of Xt, it is enough to prove that Xt does not vanish on 


each Ck- Using Eqs. (2.4) and (2.5) we can evaluate the components Pt and Qt 


of the vector field on each circle C^ = {/^ = 0}, thus obtaining 


( 2 . 6 ) XT\ck = rkBlc^fJ-klck 


Ck 


dx 


dfk 

dx 


Ck 


dy = TkXLRlCk 
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Since the gradient of fk only vanishes at the point pk , and the functions B and pu 
do not vanish on Ck (because all the circles in F are disjoint and no point Pj is 
contained in any circle Ck), we infer that Xt has no zeros on Ck, which is then a 
periodic orbit. 

Let us assume that the periodic orbit Ck is not a limit cycle. Since Xt is 
polynomial, it then follows that Ck must belong to a period annulus, i.e. it is 
surrounded by a family of periodic orbits. Consider a periodic orbit 7 ^ close enough 
to Ck so that it is disjoint from the set F and all the points pj. In particular, we 
have that Vrl'yt. does not vanish. Defining the 1-form ujt ■= —Qxdx -\- Pxdy, we 
have that 


(2.7) 



because 7 fc is a periodic orbit of Xt- On the other hand, using the definitions of 
Pt, Qt and Vt, we can also write 







d (In At) 


zbdTTSfc 7^ 0 , 


d{\nB) + d{\nC) 


where we have used that the functions In At and In B are smooth on 7 /j because At 
and B do not vanish there, and hence d (In A-r) = 0 = d (InS). The number 
Sk in the last equality is the number of primary cycles contained in Dc,., which 
is at least 1 (it is 1 if and only if Ck is primary). The sign of the last expression 


depends on how we orient the circle Ck- Since this formula contradicts Eq. (2.7), 


we conclude that all the periodic orbits Ck are limit cycles of Xt- Moreover, since 
the vanishing order of the inverse integrating factor Vt on each cycle Ck is 1 , it 
follows that all the limit cycles are hyperbolic mm- 

To conclude the proof of the theorem, let us show that each constant Tk can be 
chosen in such a way that the period of the limit cycle Ck is Tk- Indeed, since 


^rlck = '’'kXkRlck by Eq. (2.6), denoting by the period of the limit cycle Ck 
for the vector field Xlr., we have that the orbit Ck has period Tk if we choose 


Tk = 


Tk 


The theorem then follows. 


□ 


3. A PLANAR POLYNOMIAL VECTOR FIELD WITH PRESCRIBED LIMIT CYCLES, 

MULTIPLICITIES AND PERIODS 

We first prove that, for any configuration of cycles F, we can construct a planar 
polynomial vector field that realizes it as a set of limit cycles with prescribed multi¬ 
plicities. As in Sectionj^ the method of the proof is modeled upon the construction 
of Llibre and Rodriguez Ha- 

Theorem 3.1. Let F = {Cl, ..., C„} be a configuration of cycles and (toi, ..., m„} 
a set of positive integers, then F is realized by a planar polynomial vector field X^ 
with deg(A,^) < 2 (r -|- ’^mk), where each periodic orbit Ck has multiplicity mk- 
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Here r is the number of primary cycles in F. Moreover, Xm admits a polyno¬ 
mial inverse integrating factor, it is Darboux integrable and all its limit cycles are 
algebraic. 

Proof. Arguing as in the proof of Theorem |2.1[ we can assume that F is a set 
of n circles Ck = {fk{x,y) = 0}, see Eq. (13 for the definition of fk, where the 
first r ones are the primary cycles and no point pk lies on any Cj. 

We consider the auxiliary functions A,B,C defined in Section as well as the 
functions: 


Am := n f> 

k^l 

D := ABC, 


mfc 

k 


:= iDexp [ j , 

H := InD, 


k=l 




where A is a constant defined as 


A := ^(mfe - 1), 




which vanishes if and only if = 1 for every A: = 1, • • • , n, and 


fk 


l-TTlfc 


if LUfe > 2 , 


hk := >( (1 - ruk) 

Jnfk if TUfc = 1. 

We claim that the vector field Xm in the statement of the theorem can be defined 
using these functions as: 


Xm ■— Pm^x Qm^y 


with 


Pm ■■= -AmB^ - BF , 
dy 

Qm '■= AmB— -h BG . 

ox 


Observe that if = 1 for all k, then = A and F = 0 = G, thus we 
have Xm,\mk=i = where is the vector held constructed by Llibre and 

Rodriguez in m- 

To prove that Xm satisHes the desired properties, let us show that it is a polyno¬ 
mial vector held with deg(Am) < 2(r + that I? is a Darboux hrst integral 

and that Vm '■= AmB is an inverse integrating factor. Indeed, by direct computa¬ 
tions we obtain that Pm and Qm can be written as 


(3.1) 


Pm=PLRYlfr~" -BF, 

k^l 

n 

Qm = QLRYlfP~" +BG. 


k^l 


(3.2) 
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From these expressions it follows that Pm and Qm are polynomials of degree at 
most 2{n + '^rrik) — 1- Moreover, using the identities: 


dV _ f 1 dD G \ 9V _ f I dD F \ 

dx \Ddx^Am)' dy \D dy Am) ' 


one can easily check that I? is a Darboux first integral of Xm- It is also straight¬ 
forward to show that the vector held ^ is divergence-free in \ thus 

implying that Vm is an inverse integrating factor of Xm- 

Arguing as in the proof of Theorem |2.1[ we conclude that the circles {Ci, • • • , C„} 
are invariant under the how of Xm-, and that if Xm has a limit cycle, it has to be 
precisely one of these circles. Let us now prove that all of them are limit cycles. 

First we check that Xm does not vanish on each Ck- Using Eqs. ( |3.1[ ) and (3.2) 
we can write 


(3.3) 


Xm\Ck 


B|cfcAfc|cfc (A-f/™'‘ ^\ck) 



dx + 

Ck 


dfk 

dx 



Notice that the functions B, Afc and A + do not vanish on Ck because all the 

circles in F are disjoint, no point pj is contained in a circle Ck and A > 0 unless 
rrij = 1 for all j. In the case that all the limit cycles have multiplicity 1, it follows 
that A -I- = 1. Moreover, the gradient of fk only vanishes at the point pk, so 

we infer that Xm has no zeros on Ck, which is then a periodic orbit. 

Since the vector field Xm is polynomial, if Ck is not a limit cycle then it must 
belong to a period annulus. Let us assume that this is the case, and take a periodic 
orbit 7 fc close enough to Ck so that it is disjoint from the set F and all the points 
Pj- In particular, we have that the function Vm\-yk does not vanish. Defining the 
1-form LUm '-= —Qmdx -I- Pmdy, we have that 


(3.4) 


^=0 

^ m. 


because jk is a periodic orbit of Xm- Using the definitions of P„ 
proceeding as in the proof of Theorem |2.1[ we can also write 


Q m and Vm j and 


Vm 


dH 


AY^dh, 


j d (In A) -b d (Ini?) -I-d (InC) = ±47rsfc ^ 0 . 
Jjk ^ 


To obtain this formula we have used that the functions In A, In B and hj are smooth 
on because A, B and fj do not vanish there. The number Sk in the last equality 
is the number of primary cycles contained in Dck, which is at least 1 by definition. 


Since this formula contradicts Eq. (3.4), we deduce that all the periodic orbits Ck 
are limit cycles of Xm- 

To conclude the proof of the theorem, we notice that, by construction, the van¬ 
ishing order of the inverse integrating factor Vm on the limit cycle Ck is irik- Since 
the multiplicity of a limit cycle is equal to the vanishing order of the inverse inte¬ 
grating factor HIE], it follows that the multiplicity of Ck is rrik, in particular Ck is 
hyperbolic if and only if nik = 1. D 
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Combining the constructions in Theorems |2.1| and |3.1[ and arguing exactly as in 
their proofs, it is easy to check that the vector held X^m dehned as 


XTm •— “t” QTm^y j 


with 


Prm '■= —AmB ^ — BFt , 

oy 

3H 

Qxm ■= AmB ^ + BGt , 
ox 


realizes the set of cycles T with prescribed periods and multiplicities. Here the 
functions B and Ht were dehned in the proof of Theorem |2.1[ A^ in the proof of 
Theorem ixn and we set 

Fr := 




Gt ■= TfcAfc 


fc=l 


dy ’ 

dfk 

dx 


The constants tj, are chosen so that the limit cycle Ck has period Tk- Indeed, 
noticing that 


(3.5) 




Tmicfc = , 




if we denote by 7 ^'^(t) the integral curve parametrizing Gk of the Llibre-Rodriguez 
vector held realizing T, the constant Tk is chosen as 




dt 


Tk Jo [A+ 

thus implying that the period of Gk for the vector held XTm is Tfc. Here is 
the period of the integral curve 7 ^The theorem can then be stated as follows: 


Theorem 3.2. Let T = {Gi,... ,Gn} be a configuration of cycles, {mi, ...,m„} a 
set of positive integers and {Ti,..., T„} a set of positive eonstants, then T is realized 
by a planar polynomial vector field Xxm with deg{XTm) < 2(r + X)^fc); where eaeh 
periodic orbit Gk has multiplicity mk and period Tk. Here r is the number of primary 
eycles in T. Moreover, Xxm admits a polynomial inverse integrating factor, it is 
Darboux integrable and all its limit cycles are algebraic. 


4. A PLANAR POLYNOMIAL VECTOR FIELD WITH PRESCRIBED LIMIT CYCLES, 

PERIODS AND STABILITIES 

The goal of this section is to prove that, for any conhguration of cycles T, we 
can construct a planar polynomial vector held that realizes it as a set of hyperbolic 
limit cycles with prescribed stabilities and periods. To this end, we hrst show that 
the stability of each limit cycle Gk of the polynomial vector held Xt constructed in 
Section can be characterized in terms of the relative position of Gk with respect 
to the other cycles. In what follows, we will say that the limit cycle Gk has stability 
— 1 if it is stable and 1 if it is unstable. Since all the limit cycles considered in this 
section are hyperbolic, the interior and exterior stabilities are the same; the case of 
semistable limit cycles will be addressed in Section 

Lemma 4.1. The limit cycle Gk of the vector field Xt introduced in Seetion{^has 
stability (—1)'^'“, where Nk := card{Gj : Gk C Dcj,j k}. 
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Proof. Denoting by 7 fe(t), < € K, the integral curve of Xt whose image is the 
limit cycle Ck, we first compute the derivative of the angular variable 6j defined in 
Eq. (2.21 for each j: 


dOjilkit)) 


dt 


= XThk{t))-X0jMt)) 


(4.1) 


= 2rfeB(7fc(t))^fc(7fc(t)) 


{x - Xk)ix - Xj) + {y- yk){y - %) 


ik{t) 


{x - X^Y + {y - y^Y 

which is negative (positive) if the orientation of C/j induced by the integral curve 
7 fc(<) is clockwise (counterclockwise). The constants and the functions B and 
/ifc were defined in the proof of Theorem 2.1 Assuming that the point (xj,yj) is 


contained in the interior of the compact set > h easily follows that 
(x - xk){x - Xj) + {y- yk){y - yj) 


(4.2) 


(x - xj)"^ + {y - yY"^ 


ik(t) 


> 0 . 


Accordingly, since B is always positive over the sign of is given by 

the sign of fik over Ck- Noticing that the function fj = (x — Xj)'^ + {y — yj)"^ — 
is positive over Ck if and only if Ck is not contained in Dc^, the definition of yk 
implies that its sign is precisely (—1)'^'“. 

The stability of the limit cycle Ck is given by the sign of the following integral [2] 

fTk 

C := divXT{'^kit))dt. 

Jo 

Using the identity Xt ■ VVr = Ur divXr we obtain 


C = 

+ 



XT{lk{t)) ■ X \n A{-^k{t))dt 



XT{lk{t)) ■ VlnB(7fc(f))dt 


pTk 

Jo 


XTYlkY)) • Vln/fe(7fc(t))dt, 

where to pass to the second equality we have used that the functions In B and In fj 
with j Y fk are smooth on Ck, and therefore the corresponding integrals vanish. 
Using the definition of Xt, the identities (2.31 and the value of Xt on Ck computed 
in (2.61, after a few straightforward computations we can write 


-1 

rTk 

£ = — 

/ 

Tk 

Jo 

-1 

pTk 

+ - 

/ 

Tk . 

/o 


XT{lk{t)) • VlnC'(7fe(t)) 
XT(7fc(t))-Vln {BY[f;^)hk{t)). 


As before, the second integral in this expression vanishes because the function 
BYlj^k is smooth on Ck- Finally, from the definition of the function C we 
obtain 

^ ^ A d9Ylk{t)) ^ (-l)^'^47rsfc 

Tk “ Jo dt Tk 

where Sk G {1,... ,r} is the number of primary cycles contained in Dc^, and we 
have used the sign (—1)^'“ computed before (observe that there is a non vanishing 
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contribution to the integral above if and only if 9j is the angle whose center (xj , yj ) 
is contained in so that Eq. (4.2) holds). We then conclude that the stability 
of Cfc is as we wanted to show. □ 


This lemma proves that the stability of the limit cycle Ck of Xt is fixed by the 
configuration of cycles that we want to realize. In the following theorem, which 
is the main result of this section, we show how to modify the vector field in 
order to prescribe the stabilities of its limit cycles. The idea is to add additional 
cycles to the conhguration T to obtain a new configuration T so that Nk for the 
new conhguration has the desired sign. Since we do not want to realize the extra 
cycles r\r, we can remove them by adding a singular (zero) point over each extra 
limit cycle of the vector held Xt realizing T. 

Theorem 4.2. Let T = {C'l, ..., C„} be a configuration of cycles, {Ti,..., T„} a 
set of positive constants and {ui,..., u„} a set o/±l. Then T is realized by a planar 
polynomial vector field Xts with deg(XTs) < 2(3n + r), where each periodic orbit 
Ck is hyperbolic, has period Tk and stability Vk- Here r is the number of primary 
cycles in T. Moreover, Xts admits a polynomial inverse integrating factor, it is 
Darboux integrable and all its limit cycles are algebraic. 


Proof. As in previous sections, we can assume that T consists of circles. Take 
a circle Cn+k centered at pk of radius rk — evk ■ Recall that rk is the radius of the 
circle Ck and pk is its center. It is clear that we can take e > 0 small enough such 
that all the circles are disjoint and no pj lies on any Ck and Cn+k- We denote the 
whole conhguration by T := {C'l, • • • , C2n}- Observe that the number of primary 
cycles and their centers remain unchanged (a primary cycle Ck in T with > 0 is 
no longer a primary cycle in T, instead Cn+k will be primary, but it has the same 
center), and therefore the function B will be the same for T and T. 

Now we construct a vector held Xt as in Theorem 


2.1 


realizing the 2n cycles of 
r where each limit cycle Ck has an associated constant Tk that will be hxed later 


(cf. the proof of Theorem 2.1 for the dehnition of such a constant) and Tn+k = 1 
so that the limit cycle Cn+k has period for fc S (1,..., n}. It is obvious 


4.1 


from the dehnition of T that the parity of the number Nk dehned in Lemma 
only depends on the relative positions of Cn+k and Ck, and that (—1)^'“ = i/k for 
k = 1,... ,n. The lemma then implies that the cycles Ck have the desired stability. 
In order to remove the additional cycles Cn+k, we consider the functions: 

y) := [x - UkY + {y- bkf 

n 

Lts ■= Ik 

k=l 


where each qk := {ak,bk) G Cn+k is a point at the extra cycle Cn+k- Then the 
vector held: 

Xts ■= LtsXIt 

satishes the statements of the theorem. Indeed, since the factor Lts is positive 
on r, this set is realized by Xts as algebraic limit cycles, while the cycles Cn+k 
contain a singular point of Xts, thus becoming homoclinic connections. The factor 
Lts does not change the stability of each cycle Ck of Xt, which is Moreover, if 
Jk^it) is the integral curve parametrizing Ck of the Llibre-Rodriguez vector held 
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Xlpi that realizes f, see Sectionj^for the definition, the constant Tk must be chosen 
such that 

rriLR 

1 dt 

Jo ZLlFW ’ 

which implies that the period of Ck for the vector field Xts is Tfe- Here is 
the period of the integral curve 7 fc^(t). Finally, observe that the degree of Xts 
is 2(3n + r) because the vector field Xt has degree 2(2n + r) and the factor Lxe 
has degree 2n. Additionally, Vts '■= ABLts is an inverse integrating factor of Xts: 
and the Darboux first integral Ht of Xt is a first integral of Xts as well. This 
completes the proof of the theorem. □ 


5. The main theorem 


In this section we construct a vector field X that realizes a configuration of 
cycles r with prescribed periods, stabilities and multiplicities, thus establishing the 
main theorem of the paper. As in Section]^ we shall use Vk € {—1,1} to denote 
the interior stability that we want to prescribe for the limit cycle Ck (negative 
means stable in the interior and positive is unstable). Observe that the exterior 
stability of Ck is determined by Vk and the multiplicity as (—In 
the following lemma we show that the interior stability of each limit cycle Ck of 
the polynomial vector field Xttu constructed in Section can be characterized in 
terms of the relative position of Ck with respect to the other limit cycles and the 
set of multiplicities {toi, ... ,to„}. In the case that irik = 1 for all k we recover 
Lemma 14.11 

Lemma 5.1. The limit cycle Ck of the vector field XTm constructed in Section^ 
has interior stability Vk = (—1)™''+'^'=+!, where 

Mk ■■= ^ ruj 
J&Mk 

and Mk ■= {j € {1,. .. ,n} / Ck C Dc^J fc}. 


Proof. Consider a closed curve Ck in the region bounded by Ck and close enough 
to it so that Ck is disjoint from all Cj and pj. Since Ck is a limit cycle of XTm we 
can assume that Ck is transverse to the integral curves of XTm at each point. The 
interior stability of Ck is determined by the sign of the flux through Ck of XTm, 
which is given by 


(5.1) 


Fluxk '■= 



XTm ■ nds , 


where n is the unit normal vector on Ck pointing outwards and s parametrizes the 
curve Ck in the positive direction (i.e. counterclockwise). 

The sign of the flux (5.1) can be easily computed using the 1-form uiTm ■= 
—Qrmdx + PTmdy and the inverse integrating factor VTm '■= AmB of XTm- Indeed, 
noticing that Vrm does not vanish at any point of Ck, and observing that 


WTm 

Vrm 


-dHT-^Yl 

i=i 
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see the definitions of all the involved functions in Sections and it is obvious 
that the sign of FluXk is the same as the sign of 

FluXk ■■= VTm{Po) [ = - VTmiPo) [ dHx = dTTSfcVrmbo) , 

JCk JCk 

where po is any fixed point on Ci~. The second equality follows from the fact that 
each term does not contribute to the integral because fj does not vanish at any 

point of Ck- For the last equality we have used the expression of J dFlx obtained 
in the proof of Theorem 


2.1 


taking into account that the curve is positively 

oriented. 

Since the interior stability of Ck is given by — sign{Fluxk), the formula above 
implies that Vk = — sign(Vrm(G)) = — sign(Am(C'fc)). Then, arguing as in the 
proof of Lemma |4.1[ but taking into account the multiplicity, we easily obtain the 
desired expression for Vk- □ 


In the following theorem we show how to modify the vector field Xx-m in order 
to prescribe the interior stabilities of its limit cycles. The idea is the same as in the 
proof of Theorem |4.2[ we add additional cycles to the configuration T to obtain a 
new configuration T so that Pk for the new configuration has the desired sign. To 
remove the extra cycles r\r, we add a singular point over each extra limit cycle of 
the vector field Xxm realizing T. 

Theorem 5.2. Let T = {Ci,..., C„} be a configuration of cycles, {mi ,..., m„} a 
set of positive integers, {Ti, .. ., T„} a set of positive constants and {vi, ..., z/„} a 
set of±l. Then T is realized by a planar polynomial vector field X with 

deg(X) < 2 I 2{N - n) + r + y^ mk 

\ k=l 

where each periodic orbit Ck has multiplicity mk, period Tk and interior stability 
Vk- As usual, r is the number of primary cycles in T, and N := n + ni + 2n2 where 
ni := card{k : mk is odd} and n 2 ■= card{k : mk is even and Mk is odd}. 
Moreover, X admits a polynomial inverse integrating factor, it is Darboux integrable 
and all its limit cycles are algebraic. 


Proof. As usual, we assume that T consists of circles. Let us define a new 
configuration of circles L := L U [Cn+i, ■ ■ ■, Cn} following these rules: 


• If mk is odd, we add a concentric circle of radius rk — VkS. 

• If mk is even and i/fc = — 1, we add nothing. 

• If mk is even and = I, we add two concentric circles of radii ± e. 


Here e is a small enough constant so that all the circles in f are disjoint, and disjoint 
from pk. Using Lemma [5T] it is easy to check that N = n + ni + 2n2. 

Now, we construct a vector field XTm as in Section that realizes the con¬ 
figuration r, with TOfc, k G {1,. •. ,n}, the multiplicity we want to prescribe and 
m„_i_j = 1, j S {1 ,... ,N — n}. Moreover, the constant Tk corresponding to each 
limit cycle Ck will be fixed later, while we take = 1. A simple argument 

using Lemma 5.1 


implies that the interior stability of each limit cycle Ck of Xxm 


is precisely Vk- 
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Taking an arbitrary point qj := {aj,bj) G Cn+j for every extra circle, we define 
the function L as in the proof of Theorem \4.2\ i.e. 

N-n 

L := Ij , 
i=i 


where Ij := {x — ajY + (y — and the vector field 

X := LX-Tm ■ 


Since the factor L is positive on T, this set is realized by X as algebraic limit 
cycles, while the remaining cycles Cn+j C T contain a singular point, thus becoming 
homoclinic connections. The factor L does not change the interior stability of each 
cycle Cfe, which is then i>k by construction, nor the multiplicity nik- The constants 
can also be chosen such that each limit cycle Ck of X has period Tk. More 
precisely, since the vector field Xxm on each limit cycle Ck can be written as in 
Eq. (3.5), we conclude that 


Tk = 


1 

n 


rpLR 

C^k 


dt 


tN 






where we are using the notation introduced in Section 

Finally, an easy computation shows that X is a polynomial vector field of degree 
as in the statement of the theorem, and V := AmBL is an inverse integrating 
factor. Using the functions defined in Sections and it is ready to check that 
the function 

N 

AtBC exp A Tj hj ^ 

is a Darboux first integral of the vector held X. This completes the proof of the 
theorem. □ 


Remark 5.3. We observe that the vector field X in the proof of Theorem \5.^ can 
he constructed without including the functions fn+j, j £ — n}, in the 

quantities Ht, Ft an Gt appearing in the definition of XTm- 
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